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Abstract 

'Hypergeometric states', which are a one-parameter generaHzation of binomial 
states of the single-mode quantized radiation field, are introduced and their non- 
classical properties are investigated. Their limits to the binomial states and to the 
coherent and number states are studied. The ladder operator formulation of the 
hypergeometric states is found and the algebra involved turns out to be a one- 
parameter deformation of su{2) algebra. These states exhibit highly nonclassical 
properties, like sub-Poissonian character, antibunching and squeezing effects. The 
quasiprobability distributions in phase space, namely the Q and the Wigner func- 
tions are studied in detail. These remarkable properties seem to suggest that the 
hypergeometric states deserve further attention from theoretical and applicational 
sides of quantum optics. 
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1 Introduction 



The number and the coherent states of the quantized radiation field play important roles 
in quantum optics and are extensively studied [0]. The binomial states (BS, also called 
intermediate number- coherent states) introduced by Stoler, Saleh and Teich in 1985 
interpolate between the most nonclassical immber states and the most classical coheTent 
states, and reduce to them in two different limits. Some of their properties P, ^ 3, 
methods of generation ||^, ^, as well as their interaction with atoms 0, have been 
investigated in the literature. The BS is defined as a linear superposition of number 
states (|n), n = 0, 1, . . .) in an M + 1-dimensional subspace 

M 

l^,M) = E/3f(r^)|n), (1.1) 

n=0 

where 1] is a. real parameter satisfying < r/ < 1 (probability), and 



M 

n 



nl/2 



:i.2) 



The name 'binomial state' comes from the fact that their photon distribution | {n\ri, M) p 
= \l3n{'n)\^ is simply the binomial distribution with probability rj. In the two limits 
?7 — » 1 and ?7 — > it reduces to number states: 

|1, M) = |M), |0, M) = |0). (1.3) 

In a different limit of M ^ oo, t] with rjM = fixed [a real constant), \ri, M) 
reduces to the coherent states (not the most general ones, only those with real amplitude 
a), which corresponds to the Poisson distribution in probability theory 0. It is well 
known that the binomial distribution tends to the Poisson distribution in the above 
limit0. The notion of BS was also generalized to the intermediate squeezed states 
and the number-phase states |^, as well as their g-deformation |jlO|. In particular, in 



a previous paper 0] we derived the ladder-operator form the BS and on this basis 
we generalized the BS to the generalized BS which possessed the number and squeezed 
states as limits in the framework of Lie algebra su{2). 

In the present paper we shall propose a one-parameter generalization of the binomial 
states, the hyper geometric states (HGS), whose photon distribution is the hypergeo- 
metric distribution in probability theory [0. It is well known that the hypergeometric 
distribution tends to the binomial distribution in certain limit. This leads to the reduc- 
tion of HGS to the BS in the same limit. Some mathematical properties, such as the 
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equivalent ladder operator definition, related algebraic structures, will be formulated. 
It is interesting that the algebraic structure is a well-investigated generally deformed 



oscillator algebra ||T2[, which reduces to the universal enveloping algebra of Lie alge- 
bra su{2), the algebraic structure characterizing the binomial states. The nonclassical 
properties of the HGS, the photon statistical properties, sub-Poissonian distribution, 
antibunching effect and the squeezing effect will be investigated in detail. Two well 
known quasi-probability distributions, the Q-function and Wigner function will be eval- 
uated to study the nonclassical properties. It will be shown that the HGS exhibit highly 
nonclassical behaviour. 

2 HGS and basic properties 
2.1 Definition 

The HGS is defined as a linear combination of number states in an (M + l)-dimensional 
subspace 

M 

\L,M,v) = j:H:f{v,L)\n), (2.1) 

n=0 

where the probability rj is a real parameter satisfying < rj < 1. Lisa real number 
satisfying 

L>msix{Mr]-\M{l~r])-^}, (2.2) 

and 



Lr^ \ ( L{l-r,) 
n j \ M-n 



L 

M 



a 



a{a — 1) ■ ■ ■ (a — + 1) fa 



(2.3) 

1. (2.4) 



Note that in Eq.( p.4| ) the real number a is not necessarily an integer. 

The name of HGS comes from the fact that the photon distribution \ {n\L, M, r/)^ 

Kn|L,M,r/)p=[iff(ry,L)p, (2.5) 

is the hypergeometric distribution in probability theory [^]. (For a background, see the 
Appendix.) We remark that in the case of M = 1, the HGS \L, l,ri) is L independent 
and is equal to the binomial state \L, l,r]) = \l,f])- 

It is well known that the hypergeometric distribution tends to the binomial distri- 
bution in the limit L ^ oo. So, correspondingly, the HGS tends to the BS in this 
limit 

|L,M,r7)^-^|M,r/). (2.6) 
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This fact can be verified directly. Furthermore, the BS go to the number states and 
coherent states, the latter corresponds to the Poisson distribution, in certain limits. So 
the HGS reduce to the number and coherent states in these limits: 



|L,M,r/)^-=^|M,r7) 



|M), when r] ^ 1, 

|0), when rj ^ 0, (2.7) 
when M —>■ oo, rj with finite rjM = a. 



It is easy to show that the HGS are normalized by using the following identity 



M 

E 

n=0 



a 
n 



/9 

M -n 



a + 13 
M 



{21 



which can be obtained by comparing the power series expansion of (1 + t)"(l + 1)^ 
(l + t)"+^. 

2.2 Ladder operator form 



In a previous paper ||rT[ we have shown that the BS admit the ladder operator formula- 
tion, namely, they are characterized by the following eigenvalue equation 



(2.9) 



where = \/M — N a is the raising operator of the Lie algebra 5^(2) via its Holstein- 
Primakoff realization. Hereafter and a are the creation and annihilation operators of 
the photon and = a^a is the number operator. So we naturally expect that the HGS 
satisfy a generalized eigenvalue equation and the algebra involved is a deformation of 
5^(2). To this end, we suppose that the HGS satisfy an eigenvalue equation 



[/(iV)+^?(iV)a]|L,M,r/) = A|L,M,r/), 



(2.10) 



in which A is the eigenvalue to be determined. Inserting (|2.1|) into ( |2.10|) and comparing 
the coefficients, we obtain 



/(M) = A, 
(A - f{n)) 
Lrj 



(2.11) 



Lr/ \ / L{l-r]) 
n ) \ M-n 

I ( Hi-v) 



n + 1 ) \ M -n-l 
From (|2.11|) and ( p.l2|) , we have 



in + l) 



g{n), (0<n<M-l). (2.12) 



/(M) - f{n) 



Lrj — n 



L{1 - r]) - M + n + 1 
4 



y/ M — ng{n). 



(2.13) 



Observe that in the right side of the above equation M and n appear in the form M — n. 
Requiring that Eq. ( |2.10| ) reduces to ( p.9| ) in the hmit L oo, we obtain 



f{N) = ^N, g{N) = ^ 



L{l-r]) - M + N +1 



VM -N. 



Lri-N 

Substituting (|2.14| ) into (|2.10|) we arrive at the ladder operator form of HGS 



(2.14) 



\L,M,r]) = ^M\L,M,r]). (2.15) 



It is easy to see that (|2.15|) reduces to ( |2.9|) in the hmit of L —>■ oo for finite M and A^. 

Here we would like to remark that the operator in the left side of (|2.15|) is an (M + 
1) X (M + 1) matrix and it generally has M + 1 eigenvalues and eigenstates. The HGS 
is only one eigenstate of the eigenvalue y/rjM. 



Let us examine the algebraic structure involved in ( |2.15| ). Define A{L,M) as an 
associative algebra with generators 



A 



V 



'L{1 -7]) - M + N + 1^ 



M 



Jt,, = {Al,)\ (2.16) 



^1 — 7] J \ Lt] — N 

Then it is easy to verify that these operators satisfy the following commutation relations 



[N,A 



+ 4± 



Ai,A-^ = FiN), Al,A\, = F{N + 1), 



(2.17) 



where the function F{N) 



FiN) 



ri[L{l -ri)-M + N]N{M - N + 1) 



(2.18) 



{1 - r]){L7] - N + 2) 

is non-negative for < < M. This algebra A{L, M) is nothing but the generally 
deformed oscillator algebra with the structure function F{N) [0. In terms of the 
generators of A{L, M), Eq.( |2.15| ) can be rewritten as 



^N + Jl-r]A 



\L,M,r]) = ^M\L,M,r]). 



(2.19) 



It is interesting that this algebra A{L, M) is an L-deformation of su{2) in the sense 
that it contracts to the universal enveloping algebra of the Lie algebra su{2) in the limit 
L ^ oo with M and finite 



A 



-|- L— »oo 
M ^ 



(2.20) 



This means that the ladder operator form reduces to that of the BS. 
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3 Nonclassical Properties 



In this section we turn to the nonclassical properties of the HGS. 

3.1 Mean photon number and fluctuation 

The mean photon number in the HGS is obtained as 

(A^) = (L, M,r]\N\L, M, r/) = Mr/. (3.1) 

It is interesting that it is independent of L and therefore it is exactly same as that of 
the BS. However, the mean value of A^^ depends on L 

(Ar2) = ^ . (3.2) 

Then the fluctuation of the photon number is 

(AAT^) ^ (iV^) - (iV)^ = r/(l - = {AN')ss^^, (3.3) 

where the {AN'^)bs = vi^ ~ ifjM is the corresponding fluctuation of the BS. We find 
that this fluctuation is always weaker than that of the BS since the factor W{L, M) = 
[L — M) /(L — 1) is always smaller than 1 except for M = 1 and the limit L — oo. Let us 
go into some detail of the factor W{L, M), which is referred to as the weakening factor. 
We shall see that this factor lies between one half and 1 

1<W{L,M)<1. (3.4) 

In fact, as a function of L for fixed M, W{L, M) is an increasing function and the 
smallest W{L, M) corresponds to the smallest L, which is 2M for r] = 0.5. In this case 
the weakening factor is rewritten as W{M) = M/{2M — 1) which is always larger than 
1/2. 

So, in comparison with the BS, the fluctuation is reduced in the HGS. For large M 
with 1] = 0.5, the fluctuation is only about half as much as the binomial states. This is 
an important feature of the HGS. 

3.2 Sub-Poissonian distribution 



Let us introduce Mandel's Q parameter [13| defined by 



{AN') - (N) 

Q = , (3.5) 
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which measures the deviation from the Poisson distribution (the coherent state, Q = 0). 
If Q < (> 0), the field is called sub(super)-Poissonian, respectively. For the HGS, it 
is easy to find that 

Q = {l-r])W{L,M)-l, (3.6) 

which is generally negative since 1 — r/ < 1 and W{L,M) < 1, namely, the field on 
the HGS is sub-Poissonian. Exceptions are the coherent state limit (L ^ oo, M — > 
oo, rj with rjM = o? finite) and the vacuum state limit [L —>■ oo, r/ ^ 0). The 
extreme case is Q = — 1 since (1 — ri)W{L, M) is always positive. In the case 77 ^ 1 and 
L = (1 — ri)^^M —>■ 00, namely, on the number states, the extreme case occurs. 

3.3 Antibunching effect 

We say the field is antibunched if the second-order correlation function g^'^\0) satisfies 
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For the HGS 



J,, {a^a^aa) _ {N^-N) 



.^^^0)^^^ is, 

M Lt] — rj 



which always satisfies the condition (pj.T]). So the HGS is antibunched. In fact, the 



occurrence of antibunching effect and sub-Poissonian are concomitant for single mode 
and time independent fields. 

In the binomial state limit L 00, ^'^^■'(O) reduces to the second-order correlation 
function g^^\0)BS of the BS 

,(2)(o)-_^,(^)(0)^, = ^. (3.9) 

Since the factor {Lrj — l)/{Lri — t]) < 1, the HGS are more strongly antibunched than 
the BS. 

3.4 Squeezing effect [1,4] 

Define the quadrature operators x (coordinate) and p (momentum) by 

1 i 

x = -^ia'' + a), p=--j={a^-a). (3.10) 

Then their variances 

{Ax') = {x') - (a;)^ {Ap') = {p') - {p)\ (3.11) 
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obey the Heisenberg's uncertainty relation 



(Ax2)(Ap2) > t_ 



(3.12) 



If one of the (Ax^) and (Ap^) is smaller than 1/2, the squeezing occurs. For convenience, 
we define the squeezing indices 



(Ax2) - 1/2 



Sr) 



M - 1/2 



(3.13) 



1/2 ' ^ 1/2 ■ 

If 5*2; < {Sp < 0), there is squeezing in the quadrature x (p). Now let us evaluate these 
indices. 

It is easy to derive that 



a"|L,M,77) 



M—n 
fc=0 L 



-1/2 



Lr]{Lr] — 1) • ■ ■ {Lrj — n + 1) 



Lr] — n \ I L(l — r]) 
k j \ M ~ n - k 



nl/2 



\k) 



(3.14) 



for n < M and M, rj) = for n > M. In particular, for n = 1, 2, we write 



M-l 



M~2 



\L,M,r^) = ^ Hk\k), a'\L,M,r^) = ^ Hk\k), 



(3.15) 



fc=0 



k=0 



where 



L 

M 



Lr,-l\ L{l-r,) 
n \ M -1-k 



, (0<A;<M-1), 



(0 < A; < M-2). 



M 



Lr]-2 \ / L{l-r]) 
n / I M -2-k 



(3.16) 



In terms of Hn and Hn, we can obtain the squeezing indices as 



M-2 



n=0 



■A/-1 



n=0 



M-2 



Sp = 2Mr]-J2 HnHn, 



(3.17) 
(3.18) 



n=0 



in which Hn = {j], L) in (pl3|). We have also suppressed the rj and L dependence in 
Hn and Hn- 

Figures 1 and 2 are plots showing how the Sx depends on the parameter L and rj. In 
each case, different values of M (5 and 50) are chosen. From these plots we find that: 
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(1). Dependence on L (Fig.l): When or 1 — is small, \Sx\ is always larger than 
that of the BS. The smaller L, the larger 15*^.1. However, the difference from the BS is 
also small. This is understandable since in this case, L must be much larger than M 
due to the condition (|2.2|), and therefore the HGS are close to the BS. When rj is around 



0.5, the L can be closest to M (two times), and the HGS are far different from the BS. 
In particular, when L is small (close to 2M), Sx changes drastically in comparison with 
those of the BS. In general the squeezing is great for large M (50, see Fig. 1(b)) and it 
decreases for small M (5, see Fig. 1(a)). 

(2) . Dependence on r] (Fig.2): We have chosen r] = 0.25 ~ 0.75 for M = 5 and 
rj = 0.05 ~ 0.95 for M = 50. Similar to the BS, the squeezing increases as r] increases 
to a maximal point, and then it decreases rapidly. This similarity is easy to understand. 
In order to have a wide range of rj, L is much larger than M. In this case, the HGS are 
closer to the BS. We find that the HGS exhibits stronger squeezing than the BS. 

(3) . Dependence on M: From Fig.l we find that when rj is around 0.5, and for small 
L, the squeezing is weaken for small M (5) than large M (50). From Fig.2 we conclude 
that large M has wider and stronger squeezing range of 1] than the small M. The HGS 
has wider squeezing range of t] than the BS. This can be seen from Fig.l with t] = 0.72 
(a) and r] = 0.923 (b). 

4 Q and Wigner functions 

The quasi-probability distributions [|I5| in the coherent state basis turn out to be useful 
measures for studying the nonclassical features of the radiation field. In this section we 
shall study the Q and Wigner functions. 

4.1 Q-function 

We start with the Q{j3) function 

Q(/?) = -|(/5|L,M,r^)P, (4.1) 

TT 



where is the coherent states of the field. Substituting the HGS into ( [4.1| ) we obtain 
the Q-function as follows 



QiP) 



TT 



M 



n=0 



(4.2) 



Here /5 is a complex c-number P = x + iy, with {x,y = p) corresponding to the two 
quadrature operators x and p in ( |3.10| ). 
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Now we would like to investigate numerically the changes of Q-function for different 
L, M and rj. Fig. 3 are plots of Q-function of HGS for different L, for fixed M = 5, r/ = 0.5. 
When L — s> oo, the HGS is in fact a binomial state and its Q-function is shown in 
Fig. 3(d) (see also [^). Then we choose finite L values 40 (Fig.3(c)), 20 (Fig.3(b)) and 
10 (Fig. 3(a)) (note that 10 is the smallest allowed value of L for M = 5,1] = 0.5). We can 
see clear deformation of the Q-function. At first sight this deformation pattern appears 
similar to that of Q function with respect to rj given in Fig. 4 of [Q. However, they are 
essentially different: Increase in rj brings the gain of the energy or the mean photon 
number as given (|3.1| ), while the changes of L does not correspond to any change of the 



mean energy due to Eg . (|3rT|) . Fig. 3(e) is the Q-function of HGS for r] = 0.9, L = 50 and 
M = 5. This Q-function is almost the same as that of BS (see Fig.4(c) in Ref. 0]), as 
expected. 

From the Q-functions we can also study the squeezing properties by examining the 
deformation of their contours. As before we pay our attention to the case rj = 0.5 and 
explain the drastic changes of squeezing for small L. Fig. 4 are the plots of contours of 
Q-function for M = 5, 1] = 0.5 and L = 10, 28 and oo (binomial states). We find that, 
when we decrease L, the contour is first squeezed (Fig.4(b)) in the x direction until a 
maximum squeezing is reached. Then the contour deforms to the shape of an ear, which 
occupies a wider range in the x direction and the squeezing is reduced (Fig.4(a)). From 
this approach we can also explain the drastic increase of squeezing for M = 50 and 
f] = 0.5. In fact, when L becomes smaller, the shape of the contour is compressed, which 
corresponds to the strong squeezing. And this change continues until the smallest value 
of L allowed for rj = 0.5, ie. L = 100. However, in contrast to the case M = 5, the case 
M = 50 does not give rise to the shape of an ear (Fig.4(d)). 

4.2 Wigner function 



The Wigner function in the series form is defined as 



r) OO 

W{P) = -J2i-lf{P,k\p\P,k), (4.3) 

k=0 

where k) = D{P)\k) = exp(/3a^ — P*a)\k), p is the density matrix (projector on the 
HGS) and takes the form 

p=\L,M,r]){L,M,r^\ (4.4) 
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for the case in hand. It is easy to compute that 



oo M 



(4.5) 



fc = 71=0 



Here Xnk{l3) = {n\D{(3)\k) are given by [T^] 



XnkiP) = < 




exp(-|/3|V2)(r)'=-"/:^"(|/3r) ifn<fc, 



(4.6) 



where £"(|/3p) are the generahzed Laguerre polynomials. 

Since the case M = 1 is simply the binomial state and its Wigner functions have been 
investigated in detail in Ref.[Q, we here consider the simplest nontrivial case: M = 2. 
Fig. 5 are some plots of the Wigner function of HGS for M = 2 and different r]: (a) 
f] = 0.2 (b) 1] = 0.5, (c) rj = 0.9 and (d) the number state |2) (77 1 and L 00). 
In each case we choose the smallest possible value of L for given rj and M to see the 
maximal contrast with the BS. We note that the case rj ^ is just the vacuum state 
and its Wigner function is simply the Gaussian centered at the origin. As 1] increases 
from 0, this Gaussian distribution continuously deforms to the Wigner function of |2), 
as shown in Fig.5. From rj = 0.2, the negative parts of the Wigner functions are very 
clearly visible and this signifies the nonclassical properties. 

Fig. 6 are two plots of the Wigner function of the BS for M = 2 and (b') rj = 0.5 
and (c') 7] = 0.9. Comparing them with those of HGS, namely, Fig.5 (b,c), we find that 
in the case 1] = 0.5 the Winger distributions of HGS and BS are markedly different: 
distribution of HGS has two negative peaks while the BS has only one. However, in the 
case of 1] = 0.9, two distributions are almost the same. This is understandable since for 
f] = 0.9, L must be very big and the HGS are very close to the BS. 

5 Conclusion 

We have shown various properties of the hypergeometric states. The relationship with 
the BS is clarified together with the coherent state and the number state limits. The 
ladder operator formulation gives an algebraic characterization of the HGS based on 
the generally deformed oscillator algebra. The salient statistical properties of the HGS 
such as the sub-Poissonian character, the anti-bunching effect and the squeezing effects 
are investigated for a wide range of the parameters. The nonclassical features of the 
HGS for certain parameter ranges are demonstrated in terms of the quasiprobability 
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distributions, the Q-function and the Wigner function. On account of these remarkable 
properties we are tempted to think that the HGS play an important role in quantum 
optics. Surely they deserve further investigation including the method of generation. 
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Appendix. Hypergeometric Distribution 

Consider a pot containing Li red and L2 black balls. A group of M balls is chosen 
randomly. Then the probability qn that the chosen group contains exactly n red balls is 
given by the hypergeometric distribution 

-1 

n l\M -nl \M ; 



It is easy to see qn = \H^{r], L)p ( |2.5| ) for rj = Li/L. The name is explained by the fact 
that the generating function of the above distribution can be expressed in terms of the 



hypergeometric functions Obviously the above distribution tends to the binomial 



distribution with probability rj in the large L limit. 
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Figure 1: Squeezing index of HGS as a function of L, for (a) M = 5 and (b) M — 50. 
The T] values are indicated in the figures. 




Figure 2: Squeezing index of HGS as a function of rj and for different M values: (a) 
M = 5 and (b) M = 50. The L values are shown in the figures. 
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Q(x, 



Q(x, 




Figure 3: Q(/?)-functions of HGS. p = x + iy. (a) t] = 0.5, L = 10, (b) r/ = 0.5, L = 20, 
(c) 7] = 0.5, L = 40, (d) 7] = 0.5, L — > cx) (binomial state). Those four figures show how 
Q-function depends on L. Case (e) corresponds to ?7 = 0.9 and L = 50. In all the cases 
M = 5. 15 




-2 -1 1 2 3 4 -2 -1 1 2 3 4 




ure 4: Contours of (5(/3)-functions of HGS for M = 5, 1] = 0.5 and (a) L = 10, (b) 
= 28 and (c) L = oo (Binomial state); (d) M = 50, r] = 0.5 and L = 100. P = x + iy. 
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Figure 5: (3 = x + iy. Wigner function W{(3) of HGS for M = 2 and different rj and L: 
(a) rj = 0.2, (b) r] = 0.5, (c) r] = 0.9 and (d) r] = 1 (the number state |2)). The smallest 
possible value of L is chosen for each rj to show the maximal contrast with BS. 
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Figure 6: Wigner function of BS for M = 2 and (b') r] = 0.5 and (c') r] = 0.9. 
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